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By extending the Orthodox theory of Coulomb blockade to include the nano-scale effects of envi- 
ronmental electric fields, we show that tunneling spectra of metallic nanoislands contain information 
not only on the energies of the electronic levels but also on their wavefunctions near the surface of the 
island. This fundamental additional information is predicted to take the form of new observable phe- 
nomena beyond the scope of standard Orthodox theory: the tunneling resonances are renormalized 
and show surprisingly strong Gaussian fluctuations, level bending and avoided crossings. 
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I. INTRODUCTION 

The Orthodox theory of Coulomb blockade^ is a cor- 
nerstone of the present understanding of charge trans- 
port in single-electron devices^. In this theory the de- 
vice (a small metal, semiconductor or fullerene island 
and electrodes) is modeled using an equivalent electric 
circuit with macroscopic components (capacitors and re- 
sistors). Electronic interactions that control the passage 
of charge through the island are accounted for classically 
by an electrostatic charging energy U c . The quantum 
nature of the island is incorporated through the discrete- 
ness of its energy levels. The statistical properties of the 
electron wavefunctions and energy levels are usually de- 
scribed within the frame of random matrix theory^^ or 
by exact diagonalization of a microscopic Hamiltonian 7 . 
The average energy level separation S and U c set the en- 
ergy scales of the transport problem. The Orthodox the- 
ory works impressively well for mesoscopic islands. Nev- 
ertheless, recently, a new generation of nanoscopic metal- 
lic islands (nanoislands) has yielded surprising tunneling 
spectroscopy results&2ii2iii which have led theoreticians 
to extend the Orthodox theory in order to successfully 
accounliiSiiiiiii^ for some of the novel observations. 

In this paper, the Orthodox theory is extended further 
in order to explore the microscopic effects of environ- 
mental electric fields on the tunneling spectroscopy of 
these nanoislands. These effects are due to penetration 
of the fields into the island, therefore they become more 
significant as the size of the island is reduced. To our 
knowledge they have not been considered previously in 
the experimental or theoretical literature. Our paradigm 
for studying them will be the transistor introduced by 
Ralph, Black and Tinkham 11 , although environmental 
fields should affect all nanoscopic single-electron tran- 
sistors. We show that environmental fields result in an 
electrostatic renormalization of the electronic states of 
the nanoisland that is beyond the scope of the standard 
Orthodox theory. This renormalization exhibits nearly- 
Gaussian fluctuations due to the stochastic nature of the 
electron wavefunctions. Furthermore, we predict the oc- 
currence of avoided crossings and level bending due to the 
coupling between electronic levels induced by the fields. 
We show that these phenomena are reflected directly in 



the behavior of tunneling resonances; this establishes a 
platform for probing the environmental coupling experi- 
mentally. The fluctuations in the slopes of the tunneling 
resonances vs. gate voltage that we predict are surpris- 
ingly strong, comparable in size to the unrenormalized 
slope predicted by Orthodox theory. We demonstrate 
that by observing these fluctuations it is possible to as- 
sess the amplitudes of the electronic wavefunctions at the 
surface of an ultra-small conductor. Moreover, the statis- 
tics of the wavefunctions at the surface of the nanopar- 
ticle can be separated from those in the bulk and mea- 
sured. No way to measure these important properties 
of nanoparticle electron wave functions has been known 
until now. We conclude that a novel fingerprint of the 
electron wavefunctions at the surfaces of nanoscopic con- 
ductors can be found in tunneling spectroscopy experi- 
ments. 



II. ENVIRONMENTAL ELECTRIC FIELDS 

We now introduce the environmental electric fields. 
The model geometry^ is shown in Fig. da): Three 
macroscopic electrodes labeled i = L and R (the source 
and drain contacts) and G (the gate), are separated 
from the central nanoisland (P) by thin insulating lay- 
ers of thickness d l , and are connected to voltage sources 
Vl i Vr, and Vq which drive the island to a potential 
V P = (Cs)- 1 [Q P - (C L - C R )V/2 + C G V G ] in the Or- 
thodox theory. Qp is the charge on the nanoisland. Cx, 
Cr, and Cq are capacitances between electrode i and the 
nanoisland, and Cs = Cl + Cr + Cq. V = Vr — Vl is 
the bias voltage with — Vl = Vr = V/2. The difference 
between V* and Vp sets an electric field at the nanois- 
land surface (fij) that faces electrode i. The electric field 
has a finite penetration depth (d s ) into the conductor (in 
the standard Orthodox model d s = 0) that causes the 
electrostatic potential within the nanoisland to deviate 
from V P . This deviation (AV* = V* - V P ) is sketched 
in Figure db) . Assuming that the electric field within 
the metal nanoisland is screened exponentialljii&, AVI 
at a distance z along the direction normal to f2j can 
be estimated as AV* = a i (Xd i + 2)" 1 (Vi - V P )e~ Xz i2. 
cti = e(\d l + 2)/(\d l + 2e) is a polarization enhancement 
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factor, where e is the dielectric constant of the insulat- 
ing material, and A = l/d s (~ lA _1 for Ali&). Due to 
the effective screening (small d s ) in metals, only atoms 
near the surface are significantly perturbed by the elec- 
trostatic potential energy U — — AV^,- However, in 
the nanoscopic regime the substantial surface-to-volumc 
ratio (iVa/iVv ~ 1) of the nanoisland warrants the inclu- 
sion of U in the calculation of its electronic structure. U 
causes what we call the microscopic electrostatic renor- 
malization of the energy levels, that is not included in 
the standard Orthodox theory. This renormalization is 
in addition to the bulk or macroscopic one due to — eVp 
that is included in the standard Orthodox theory. 

III. RESULTS AND DISCUSSION 

We begin by exploring the effects of U analytically 
within perturbation theory. We then present results of 



our non-perturbative computer simulations that apply to 
specific real systems: Al nanoislands with oxide tunnel 
barriers. 



A. Perturbation theory in It 

We define the microscopic contribution to the energy 
£*(V, V G ) of level \%) to be e«(V, V G ), where q = Q P /e 
labels the charge state of the nanoisland. Let be an 
electronic eigenstate for V G — 0. For V G ^ these states 
become coupled: U a b = ("^aWl^b) ^ ®' Now SU PP 0SC 
the coupling strength lab = U ab /[£%(0,0) - ££(0,0)] is 
small for all \ip h y ). Then, to first order in U, (setting for 
simplicity V — q — 0) £°(0, V G ) is linear in V G with slope 
(see Appendix lA"|l 
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where IV'oO)! 2 is the amplitude of \tp®} a t atomic site j 
of the nanoisland. z* is the distance of site j from the 
surface f^. Because d s is smaller than the size of the 
nanoisland the leading terms in the summations corre- 
spond to atomic sites on the surfaces fii. The surface 
of the nanoisland presents atomic scale disorder. Conse- 
quently, in the absence of magnetic fields, the electronic 
structure is described by the Gaussian orthogonal en- 
semble of random matrix theory 5 ^. Therefore, |^a(j)| 2 
is a random variable that obeys a \ 2 distribution^. Thus 
the slope fluctuates from state to stated reflecting 
the microscopic details of the electron wavefunctions at 
the surfaces of the nanoisland. Due to the form of Eq. 
||TJ the fluctuations are nearly-Gaussian according to the 
central limit theorem. These observations remain valid 
for arbitrary V and q. We conclude from this analysis 
that the microscopic effects of the environmental electric 
fields lead to a fluctuating renormalization (microscopic 
electrostatic renormalization) of the electronic levels of a 
conducting nanoisland. It should be noted that in the 
standard Orthodox model (d s = 0) S™ is zero. 

Consider now the implications of this microscopic 
effect for the tunneling spectroscopy of the nano- 



transistors. Suppose Qp — at V = and fixed V G , 
and then, after increasing the bias, an electron tunnels 
in from electrode L (L — > P) so that Qp = — e. If 
\ip a ) is being populated by the tunneling electron the 
transition is favorable only if the final state energy of 
the electron is lower than £p, the Fermi energy of elec- 
trode L: £ F > £-(V,V G ) = e-(V,V G ) + U c - [{C R + 
C G /2)/C s ]eV~(C G /C s )eV G , where U c = e 2 /2C s is the 
charging energy 1 . (Note that the macroscopic renormal- 
ization due to Qp (through U c ), V and V G constitutes the 
Orthodox model.) A similar restriction applies to tunnel- 
ing transitions from the neutral nanoisland to electrode 
R. These conditions for tunneling set the threshold bi- 
ases (V th ), the minimum values of V that satisfy them 
for a given \ijj a ). As V is swept at fixed V G , each time a 
new V th is reached a peak (a threshold resonance (TR)) 
appears in the differential conductance dl/dV—. We now 
analyze the effect of V G on the 77?, focusing on L — > P 
transitions. In the perturbative regime we find (see Ap- 
pendix [HJ 
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where J™ = 5e-(V,V G )/(e5V)SL. Equation © shows 
that the positions of the resonances in the dl/dV spec- 
trum shift linearly with Vg- The first term on the right 
in Eq. @ is the (macroscopic) slope predicted by stan- 
dard Orthodox theory^; the other terms are due to mi- 
croscopic effects. This analysis shows that the expected 
macroscopic slope of the tunneling resonances is micro- 
scopically renormalized. Furthermore, the slope of each 
resonance in the spectrum is different as a consequence of 
the fluctuations. Therefore it should be possible to study 
the coupling of the nanoisland's electronic levels to the 
environment by tunneling spectroscopy, and thus obtain 
experimental information on the electron wavefunction. 



B. Non-perturbative computer simulations 

So far our exposition, although perturbative, has been 
general. To understand the effects of the electric fields on 
the electronic structure and tunneling spectroscopy be- 
yond perturbation theory, we now consider Al nanoislands 
coated with Al-oxide within a tight-binding modeli with 
the on-site orbital energies of the isolated neutral Al/ Al- 
oxide nanoisland adapted to include U. The presence 
of the oxide coating introduce surface disorder in these 
nanoislands, as explained in Ref. (see also Ref. I22T) . We 
consider disc-shaped nanoislands of volume V ~ 17rtm 3 , 
and d R = lOA, d L = 15Aand d G = Ad L & The latter 
choice assures that tunneling from the gate is suppressed. 
The parameters for two such systems are shown in Table 
ID The Fermi energies of the neutral isolated nanoislands 
(E F ) differ due to their differing structural disorder and 
geometries. The latter are also responsible for the dif- 
fering capacitances. The charging voltage Vjj = e/Cs is 
the upwards (downwards) shift of Vp as one electron is 
removed (added) from (to) the nanoisland; the degener- 
acy gate voltage Vq = ±e/(2C G ) is, within the standard 
Orthodox theoryi, the value of Vg at which an electron 
is added (+) to or removed (— ) from the nanoisland at 
V = 0. 

Figure [21 shows results of the exact diagonalization of 
the tight binding Hamiltonian for e a (V, Vg) for a few en- 
ergy levels of transistors A and B under different bias 
and gate conditions, at q — 0. Vg ranges from approx- 
imately Vq to Vq . V lies within the first step of the 
Coulomb staircase (at Vg = 0) of the nanoisland. Let us 
first focus on the effects of Vg at V = 0. The striking 
features are: i) fluctuating renormalization of the slopes 
of the energy levels due to the applied electric field, ii) 
linear dependence of some energy levels on Vg, and iii) 
level bending and avoided crossings with increasing \Vg\- 
These features are qualitatively similar regardless the val- 
ues of q or V (not shown). We now proceed to discuss 
them, i) and ii) can be understood within the perturba- 
tive argument presented above. Regarding iii) it should 
be noted that as \Vg\ increases significant enhancement 
of jab can occur due to increased coupling between lev- 
els and level proximity induced by Vg- In this case, the 



microscopic contribution to the energy levels should be 
calculated to higher order in IA. Level bending signals 
the response of the energy levels to this change in the 
electronic coupling. As \V G \ increases some energy levels 
become nearly degenerate. Since these levels are in gen- 
eral mixed by U this results in avoided crossings in the 
spectrum. Thus level bending and avoided crossings as 
a function of gate bias should be general features of the 
spectra of these nano-transistors: Electric fields should 
effectively couple quasi-particle levels (which may even 
involve electronic excitations) that belong to the same 
Hilbert space and whose wavefunctions are such that the 
strength parameter "f ao is significant. The microscopic ef- 
fect of the source-drain voltage V on the electronic struc- 
ture is also shown in Fig. [21 (for q = Vg — 0) ; there are 
again fluctuations and avoided crossings. These features 
are similar for finite values of q and Vg (not shown). The 
renormalization as a function of V, however, is smaller 
than for Vg due to the electrode geometry. 

In conclusion, we have shown that the electronic lev- 
els are sensitive to the environmental electric fields and 
different regimes can be identified depending on the 
strength of the electric coupling between levels. The 
novel effects are predicted to occur in realistic Al/oxide 
transistors at the meV scale (see Fig. 0) that is readily 
accessible experimentally. 

We now discuss the fluctuations of the slopes using the 
calculated e a (V, Vg)- Table[]]shows the average value (<S) 
of the slopes S™, taken over n = 101 consecutive levels 
distributed symmetrically with respect to E F , and the 
Orthodox slope — Cg/Ce for comparison. S for transis- 
tor A is greater than for B by a factor of 1.34. This 
is consistent with £Iq being 1.38 times bigger in A than 
in B (see Table and, consequently, more atomic sites 
being perturbed by the gate electric field. The standard 
deviation (a) of S™ 1 is roughly equal (~ 2.3 x 10~ 3 ) for 
the two transistors. Figure [31 shows the histograms for 
5™ in units of — Cg/Ce- The microscopic slopes of the 
energy levels are clearly significant fractions of the macro- 
scopic slope arising from the Orthodox model. Indeed, 
for B the microscopic and the macroscopic effects are 
roughly equal. Thus the microscopic effects of environ- 
mental fields on the electronic structure and tunneling 
resonances that we have derived here can be a very im- 
portant correction to macroscopic effects predicted by 
standard Orthodox theory. Both distributions in Fig. 
13 are nearly Gaussian (the smooth lines in Fig. [3| are 
nAs(aV2%] r x exp [-(S™ -S) 2 o-- 2 /2], As = 2 x 1CT 3 
is the bin size), as predicted from Eq. (JTJ. The inset of 
Fig. EJshows J™ 21 in units of J = -(C H + C G /2)/C E . 
Clearly, this effect is negligible when compared with S™. 
Therefore, experimental measurements of SV /SVg can 
accurately determine <S™ and hence the microscopic ef- 
fects of the environmental fields on the nanoisland energy 
levels; see Eq. |J2}. 
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IV. SUMMARY 

In summary, we have shown that coupling to envi- 
ronmental electric fields gives rise to a unique finger- 
print of the electron wavefunction at the surface of a 
nanoscopic conductor that can be observed in tunnel- 
ing experiments. While we illustrated this by consider- 
ing Al/Al-oxide nanoislands our findings rely only on the 
ability of ultrasmall conducting islands to screen electric 
fields, and their large surface-to- volume ratio. There- 
fore, other metallic nanoislands (noble and magnetic) and 
more exotic nanoscopic conductors such as metallic car- 



bon nanotubes^ may also reveal the fingerprints of their 
electron wavef unctions. We hope our results will moti- 
vate experimentalists to look for these fingerprints. 
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APPENDIX A: DERIVATION OF EQ. © 

In Section |n] we showed that the applied bias (V) and gate (Vg) voltages drive the nanoisland to a potential Vp 
that depends on the charge state (Qp) of the nanoisland, and the values of V and Vg- Applying only Vg to a neutral 
nanoisland results in Vp = (C g /Cy,)V g . In this case (V = Qp = q = 0), the electrostatic potential energy U is given 
by: 



U 



eV G {C G /Cv) -a L (Xd L + 2)~"e^ - a R (Xd H + 2)- 1 e'^ + [C^/C G - l)a G (\d + 2)'^-^ 



(Al) 



Note that Eq. IjAlfl becomes different for other Qp, V, and Vg values. 

Within first order perturbation theory in U, the microscopic contribution to the energy £°(0, Vg) of level \ip a ), 
£°(0,y G ),is: 



s° a (0,V G )= ei(o,o) + mum 

= -eV G (C G /Cx) Ej l^(j)| 2 [-a L (Xd L + 2)- 1 e~ Xz > - a R (Xd R + 2)^e~ Xz f 
+ (Cz/C G - l)a G (Xd G + 2)- l e~ Xz i 



(A2) 



where \ip®) is an electronic eigenstate for V G — 0, as mentioned in the text (Sec. IIII Af) . and the summation is 
performed over all atomic sites. 

Equation in the text is straight forwardly derived from Eq. (|A2|) . 



APPENDIX B: DERIVATION OF EQ. © 

As discussed in the text (Sec. IIII Af) . V^^ p (Vg) is the minimum (threshold) value of the bias voltage — in the 
presence of an applied gate voltage Vg — at which an electron tunnels from electrode L into an empty state \ip a ) in 
the nanoisland (P). This threshold value is given byi 

£f = e-(V[t P (V G ), V G ) + U c - [{C r + C G /2)/C^}eVlt P (V G ) - (C G /Cj;)eV G . (Bl) 
After a simple manipulation, Eq. (|B1|) leads to: 



SVftpiVc) 
SV G 



Cr 



Cr + C G /2 



i + {-Cc/c^r 1 



eg(Vltp(VG + SV G ), Vg + 6V G ) - £-(^l P (Vb), Vg) 
eSV G 



(B2) 



Then, we calculate within first order perturbation theory in IA (setting Qp = — e and taking V and Vg as variables) 
the values of e~(V, Vg)- The latter leads to the following identity: 
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s-AvttAVo + y G + sv G) - e- A vttAv G ), vo) _ Va)/{e5VG) = ST + JT mtAVc)\ {m) 



e5V G ay ' w/v a V <$Vc 

where J™ 1 = fe~ (V, Va)/(eSV). By combining Eqs. (|B2|) and (|B3|> we get: 

5VtL l(} VG) = ^ /9 [1 + S?{-C G /Cv)- l ]{l + J™[-(Cr + Cg/2)/^)- 1 ]}- 1 (B4) 
oVg C fl + C G /2 

Finally, by expanding the denominator of Eq. l|B4j) in powers of J™[— (Cr + Cq/2)/Cy] we obtain Eq. @. 
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label 


D (tim) 


h (nm) 


Nn/N v 


E% (eV) 


C R (aF) 


C L (aF) 


Co (aF) 


e/C E (mV) 


e/2C G (mV) 


5 


— Cg/Ct, 


A 


3.646 


1.62 


0.774 


8.289 


1.480 


1.057 


0.164 


59.32 


488.47 


-0.0214 


-0.0607 


B 


5.266 


0.81 


1.045 


8.281 


3.088 


2.206 


0.109 


29.65 


734.95 


-0.0159 


-0.0202 



TABLE I: Transistor parameters: Diameter (D), height (h), surface-to- volume ratio (Nn/Nv), Fermi energy of the neutral iso- 
lated grain (E F ), capacitances, and charging and degeneracy voltages. S is the average slope due to microscopic renormalization 
and —G'g/Cs is the Orthodox model slope (see text). 



FIG. 1: (a) Schematic of the model transistors. The island (black and white circles) is disc-shaped and is surrounded by the 
gate (G). White circles indicate where environmental electric fields are important, (b) Sketch of the deviation AV Z within the 
island. O' indicates its center. 



FIG. 2: Fingerprints of the electronic wave functions: microscopic contribution to the energy levels of the transistors vs. Va 
(at V = 0) and V (at V a = 0) with Q P = 0. Voltages are in mV. 



FIG. 3: Histograms (broken lines) of in units of the standard Orthodox theory slope: —Cg/Cs. Gaussian distributions 
(smooth lines) corresponding to the calculated S and a (see text). Inset: Histograms of J™ in units of Jo = —{Cr + Cg/2)/Cs 
for A (bold line) and B (thin line). 
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